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completeness and orthogonality relations is strictly proven. Different Green 
functions in the form of sums over the exact solutions are constructed. The 
Fock-Schwinger proper time integral representations of these Green functions 
are found. As physical applications we consider the calculations of different 
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I. INTRODUCTION 

Exact solutions of the relativistic wave equations (Klein-Gordon or Dirac equation) in 
an external electromagnetic field as well as Green functions of those equations are very 
important in QED with such a field. Special complete sets of the exact solutions and special 
kinds of the Green functions allow one to calculate different quantum effects, for example, 
particles scattering, pair creation and so on, in zero order with respect to the radiative 
corrections, but taking the interaction with the external field into account exactly 
They may serve also as a basis for the perturbation theory with respect to the radiative 
interaction, in which the external field is taken into account exactly (so called generalized 
Furry picture) In four- dimensional space-time the above exact solutions were studied 

for different configurations of the external field in numerous papers and books (see, for 
example, 0,0. Among all the configurations of the external electromagnetic fields, which 
admit the exact solutions in d = 4, one is especially important due to the fact that it 
corresponds to a wide class of real physical situations. It is a combination of constant 
uniform electric and magnetic fields and plane wave field. The exact solutions can be found 
if there exist a reference frame in which the constant electric and magnetic fields, and the 
wave 3- vector are collinear. First, the exact solutions of Klein-Gordon and Dirac equation 
in such a configuration were found in 0. Some of the Green functions of scalar field in the 
same configuration were studied in the papers p|. Using those exact solutions and Green 
functions different quantum effects were calculated 0. Calculation of Green functions of 
spinor field was not present in the literature. Moreover, lately, field theoretical models in 
dimensions different from d = 4 attract attention due to various reasons. One can mention 
here models in 2 + 1 dimensions which probably describe planar physical phenomenon, and 
models in > 4 dimensions in connection with Kaluza-Klein ideas (see e.g. [|1^]. That is 

2 



why in the present article we are going to study exact solutions and Green functions (both 
scalar and spinor case) in ci- dimensional QED with an external electromagnetic field, which 
may be considered as a generalization of the above mentioned special configuration in four 
dimensions. At the same time we also present new results in four dimensions. A general 
definition of such a configuration valid for any dimensions of space-time is the following: 
an external electromagnetic field with constant and homogeneous field invariants. Also, an 
important motivation to consider exact solutions and Green functions in such a configuration 
of the external electromagnetic field is the fact that formally problems in some configurations 
of an external gravitational field lead to the same kind of equations, so that the former exact 
solutions can be used in the latter case after some simple identifications. Such a reduction 
may be done, for example, for De Sitter and FRW metrics |]lT],|12], for four fermionic models 
in an external field and so on. 

In a sense, the paper can be considered as a continuation and generalization of the one 
T3| , where only an electric-like external field was considered. The latter paper can help 



the reader to better understand the logic of the exact solutions constructed in this more 
complicated present case. Thus, here we sometimes omit technical details and explanations 
which can be extracted from the above mentioned previous paper. 

The paper is organized as follows: In Sec. II we describe first the external electromagnetic 
field under consideration. Then, we present the so called in and out complete sets of exact 
solutions of Dirac and Klein-Gordon equations in the external field. The cases of even and 
odd dimensions have to be considered separately, they are essentially different. We attract 
attention to the asymmetry of the quasienergy spectrum, which appears in odd dimensions. 
We spend enough time to prove strictly the in and out classification of the exact solutions as 
well as to prove the completeness and orthogonality relations. Then in Sec. Ill we construct 
different Green functions in the form of sums over the exact solutions and present them in 
the form of contour integrals over the Fock-Schwinger proper time. Among them are the 
causal Green function, in — in and out — out Green functions, the commutator function and 
so on. One ought to say that the spinor case is treated first, even in d = 4 dimensions. As 



physical applications we consider in Sec. IV the calculations of different quantum effects 
related to the vacuum instability in the external field. For example, we present mean values 
of particles created from the vacuum, the probability for the vacuum remaining a vacuum, 
the effective action, the expectation values of the current and energy-momentum tensor. We 
calculate the latter quantities by means of both in and out sets of the exact solutions. 



II. COMPLETE SETS OF EXACT SOLUTIONS 

The Dirac equation in an external electromagnetic field with potentials A^{x) in d di- 
mensions has the form [h = c = 1) 

{V^-i^ - m) 4j{x) = , V^ = id^- qA^{x) , (1) 

where ip{x) is a 2t^l-component column, are 7-matrices in d dimensions ||14|| , 

[7^71+ = 2r7'^^ r7^'^ = diag(l,-l,-l,...), d = D + I, 



d 

^0 



and X = (x^) = (x , x), x = (x*), /i = 0, 1, . . . , D, i = 1, . . . , D. The time-independent 
scalar product of the solutions of the Dirac equation may be chosen in the conventional form 



(^,^') = y ^(x)7V(a;)c/x. (2) 

As usual, it is convenient to present iplx) in the form 

i^{x) = (V^Y + m) (pix) . (3) 

Then the functions have to obey the squared Dirac equation in d dimensions, 

(r^ -rn''- IcT^'y",.) m = , J^^, = d.AAx) - d,A^{x) , = '-[Y, Y] ■ (4) 

To construct the above mentioned generalized Furry picture in QED with an external 
field one has to find special sets of exact solutions of Eq. (p, namely, two complete and 
orthonormal sets of solution: |±'?A{n}(a^)} which describes particles (+) and antiparticles (— ) 
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in the initial time instant (x° —oo), and which describes particles (+) and 

antiparticles (— ) in the final time instant +00). According to the general approach 

such solutions obey the following asymptotic conditions 

Ho.pXx^) cV^{„}(x) = cV^{„}(x), , sgn = (, x° ^ -00 , 

Ho.pix^) ''^{z}(a;) = '^e ^i){i}{x), sgn = (, x° ^ +00 , (5) 

where (, {n} and (, {1} are complete sets of quantum numbers which characterize solutions 
^ilj{n}{x) and respectively, Ho,p, = 7°(m — Y'Pi) is a one-particle Dirac Hamiltonian 

in convenient external field gauge Aq{x) = 0; ~^e, +e are particle quasi-energies and |^£:| and 
|_£:| are antiparticles quasi-energies. All the information about the processes of particles 
scattering and creation by an external field (in zeroth order with respect to the radiative 
corrections) can be extracted from the decomposition coefficients (matrices) G 



«^(x) = +^(x)G (+|«) + _^(x)G (^1^) . (6) 
The matrices G (cl'' ) obey the following relations, 

G(cr)G(cr)VG(cr)G(cr)^ = i, 

G(+r)G(_r)VG(+|-)G(_|-)' = 0, (7) 

where I is the identity matrix. 

Let us describe the external electromagnetic field in which we are going to construct the 
exact solutions. Such a field is a generalization of the corresponding field in c? = 4, which 
is a combination of the constant uniform field with the plane wave field and which admits 
there exact solutions. One can describe such a field in c? = 4 in arbitrary reference frame 
saying that both its field invariants do not depend on space-time coordinates. In > 4 
there exist, generally speaking, more independent field invariants. Besides Ji = l/2J-'^i^J-'^'^ , 
those are all possible independent invariant combinations which may be constructed from 
the field tensor jF^,y and Levi- Civita tensor. It is easy to see that there exist [d/2] such 
invariants. We define the external field under consideration in arbitrary dimensions in the 



same manner, all its field invariants have to be constant and uniform. One can see that such 
a field is a combination of a constant uniform field F^j, and a plane- wave field f^y{nx) 

^fiu = F^u + f^lv{nx), (8) 



where is an real isotropic vector, n^n'^ = 0. It is an eigenvector of the tensor F^^, F^j^yTf = 
£n^, , and f^u{nx) is a transverse field with respect to n^, n^fnu{nx) = f^u{P'X)rf = 0. If 

all of the invariants are equal to zero (it is only possible for d > 2), only a plane- wave field 
configuration is possible. If some of the invariants are not equal to zero, then a constant, 
uniform component F^^,, is not zero. The field is free. Such a field is of special interest 
due to the fact that QED with a free classical field can be treated as exact QED (without 
external fields) but with some special (coherent) initial photon states [ITSl ,^. 

If the eigenvalue S is not zero, the external field violates the vacuum stability (creating 
particles). One can find an inertial frame where the matrix J-'^i, has a simple form 



D-l 



Fl = E [616- - 6-61. , 



fi,v {nx) = ^ {n^6^ - 71^6^^ fk (x_) , /fc(x_) 

k=l 

= 6^^ — 6^, x^ = nx = x^ — X- . 



dx- 



If d > 2 is even then 

{d-2)/2 

and if d = 2 the fields and f^^, (nx) are absent. In the case d is odd, and S 0, 

{d-3)/2 

= E H,{6l^'r' - 6^6f^-') ifd>3, and = if = 3 . 

The same formula is valued if £^ = 0, but at least one of the imaginary eigenvalues of F^^^ is 
zero. Here fk{x-) are arbitrary functions of x_. If £^ = 0, and all imaginary eigenvalues of 
F^i, are not equal to zero (all field invariants are not equal to zero) we have 



E = 0, R 



{d-l)/2 

E H,{^^-'-^^-'l /.(x_) = Oforallfc. 



In the reference frame under consideration £ = E. 

To realize the external electromagnetic field of the above form we select the following 
potentials: 



A, (x) = (x) + Al {x) + f, , 

Al {x) = Ex%^^, Af = -H, 3 = 1, 

^(x_)=Oif //^O,!^. 



(10) 



,[d/2]-l, z = l,...,D-l, 



Below we present linearly independent sets of solutions of the squared Dirac equation, 
which correspond to the particles in the initial time instance and to the antiparticles in the 
final time instant: 

</'n,r(a;±) = 0pi,ni(a:i,a;2)0p3,n2(^3,a:4)---0pd_3,n(d_2)/2(2^i?-2,a;D-i), if d is even, 

{271)^^^"^ exp{—ipD~ix^~^) if d is odd and H(^d~i)/2 = 0, 
+(pp.,nAx^ x"") = (47r)-^/2 II (qE{xl/2 - xl) - + A In (Tff_)) + 



P- — T 

qE 



T \ „ I — T 



qE 



T-^dr, 



T-^dr, 



„0 I _D 



where p_, n = (ni,n2, . . . , n[d/2]-i;pi,P3, ■ ■ ■ ,P2[(d-i)/2]-i) , ^ and r = (ri,r2, . . .,r[d/2]-i) is 
a complete set of quantum numbers. Among them p_, pj are momenta of the continuous 



spectrum, rij are integer quantum numbers, ^ = ±1 and Tj = ±1 are spin quantum numbers; 
the momentum p„ is the eigenvalue of the operator 2i^|^; if qE > is chosen then the signs 
=F assigned to the functions ^4'p_^n,s,,r{x) are matched with those of the kinetic momentum 

j9_ — qEx^ at 2;_ ±oo, and 



TT 



= if /i = 0, = x'^ if = 1, . . . , D - 1, 



d 



7r_L^ = if /i = 0, D, 7r_L^ = i-^ - qA^ (x) if /i = 1, 
["^/^i-i I 0^ is even 

i=i p^_2, ciis odd , 



\qH,\{2n, + l-rj), 



H, = 



Each function 0p2j-i,nj (^sj-i, a;2j) obeys the following equations 



{-^1.23-1 + Trlaj - ^j)(i>P2,-i,nj{x2j-l,X2j) = 0, 

(7r±2j-i -P2j-i)0p2,_i,n,(a;2i-i,a;2i) = 0. 



If Hj 7^ 0, a solution of these equations is 



D-h 



(12) 



(X2j_l,2;2j) 
1/2 



\qH, 



exp < -zp2j-ia; 



2i-i 



2 V"" ■ qHj J j [V I- 71 ^ 

where 7i„^(x) are the Hermite polynomial with integer = 0, 1, . . .. If iJj = 0, the dis- 
crete quantum numbers rij have to be replaced by the momenta p2j, and the corresponding 
function has the form 

(f>P2j-i,P2jix2j-UX2j) = (27r)-^exp{-i {p2j-lX^^~^ + P2jx'^^)} ■ 

The symbol ln(=l=7f_) means the principal branch of the logarithm, ln(=l=7f_) = ln|7f_| + 
Z7r0(±7r_), while the integration paths in the r-plane, as well as the arguments 7r„, are 
shown in FIG.|I] and FIG.^ for the functions , + J and , ^ J (xJ) re- 

spectively. 
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FIG. 1. Pass of integration in ^K{x-) and +J(a;_) integrals. 
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7r_ +CXD 



FIG. 2. Pass of integration in K(x-) and J{,X-) integrals. 

One assumes the functions fk{x-) obey the requirement that (x^) and ^J(a;_) are 
analytic functions, and in particular, the functions fk{x-) vanish at ±00 quite rapidly. 

The solution has a different form if d is odd, S = 0, and all the imaginary eigenvalues of 
Ffj^i, are not equal to zero (in this case a plane- wave field is absent). 



0n,r(-^) 0pi,ni ■■■4'pd-2,n(d-i)/2i-^D-lj ^d) 



(13) 



n,riS^ ) — C exp fi^ Iqi^nrl ) ; |=F^' 



\ 



(d-l)/2 



+ E 



31 



where c is a normalization constant. 

Here v^^r are some constant orthonormal spinors, f^^r'^^.f' = ^r,r'- Equation (^) allows one 
to subject these spinors to some supplementary conditions. 
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E^v^^^r = 0, S± = 1(1 ± 7°7^), rank S± = J(,) = 21^1-^ (14) 



Rjv^^r = rjV^^r, d>A, Rj = i7^^"S^^- (15) 

If c? = 2, 3, the independent quantum number r does not appear. 

One can verify the solutions of the Dirac equation with different C,, namely, 
i'P^'j^' + m) +0p_,„,+i,r(a;) and {V^'j^' + m) +0p_,n,-i,r(a;), or (V^'j^' + m) "0p_,n,+i,r(a;) and 
{Vfj^'-i^ + m) ~0p_,n, -1.7(2;) are linearly dependent for each sign "+" or "— ". It means, in 
fact, that the spin projections of a particle (+) and an antiparticle (— ) can take on only 
J((i) values. Thus, to construct the complete sets it is sufficient to use only the following 
solutions: 

+'4'p-,n,r{x) = {V^'-f'' + m) +0p_,n,+l,r(a;) • (16) 

In particular, if d = 2,3, there is only one spin projection, with only one spinor v+ir = 
To get the case where = one has to consider a limit in the solution _,_^/'p_,n,r(2;), 



1 

V2 



such that — (qE) ^ln(=l=7r_) and > for {+ipp_^n,r{x)} whereas < for 

Using the solutions (|T3p we may construct the Dirac spinors in a slightly different way, 
;^n,r-(a;) = {V^Y + m)l (f)n,f{x) ' (17) 

+ (pnA^) = (+(Pn,+lA^) + {<lH(^d-l)/2)^(t>n-lA^) ) > ^ = ('"l, ^^2, • • . , r(rf_3)/2 , SgU (^gi^'(rf_i)/2) ) , 

where 

Rj +0n,r-(a;) = fj ~0„,j=(x), j = 1, 2, . . . , (d - l)/2 , 

and '^(f)n,f{x) are eigenvectors of the time-independent operator — V^Vi + ^a'^^F^y. 

The solutions ^ipnA^) fro^i (p!7D describe particles (+) and 'ipnA^) describe antiparti- 
cles (— ) in any time instant. They form a complete and orthonormal set of solutions. There 
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is an interesting asymmetry in energy spectrums of the such particles and antiparticles in odd 
dimensions. If all quantum numbers rij = 0, and rj = sgn {qHj) , then 7P^?/'„ f (x) = 0, and 
the ground state of the particle has energy ^eq = m, but the energy of the ground state of 



J I- 



the antiparticle is dependent of the magnetic field and different: | eo\ = \Jm? + 2min \qH 
In the d = A case the solutions similar to (|T6|) were found first in [|^ . In case E = {] they 



coincide with ones in [|16[, and if = they coincide with well-known Wolkov form [Il7| . 
In the case = 0, the solutions form a complete and orthonormal set, moreover j^'il)p_^n,T{x) 
describe particles and ~ipp_,n,r{x) describe antiparticles in any time instant. It is no longer 
if 7^ 0. However, solutions (|l^) can be used to construct new complete sets which can be 
classified by the correct way. 

One can form two complete and orthonormal sets of the solutions: (3;)} and 

^ipp-,n,r{.x) using ([16|) and additional sets 

+^p_,n,.(x) = e(7r_)(+^(a;)G(+|+)) 

-V^p_,n,.(x) = e(-7r_) (~^{x)G , (18) 

where G (^^|'''^ obey the relations (0) and, in particular, G {^\~) are decomposition coeffi- 
cients of ~ipp_,n,r{x) solutions in +ipp_,n,r{x) solutious, 

G'f+r) , , , = {+i^p.,n,r,~'ipp' ,n',r')- (19) 

\ / p^,n,r,p_,n' ,r' 

The corresponding Klein-Gordon solutions follow from (|Il|) by the replacement ug^,,(x_) = 
exp{-iln(=Fvr-)}. 

The orthonormality of all the solutions can be verified by using the following integral 
transformations: 



hoc 



±^p_,nAx) = i2nqE)-^/^ M*{pD,P-)±iJp^,nA^)dpD, (20) 

J —CO 

r+oo 

±V'pD,n,r(a;) = {2TiqEy^''^ / M{pn,P-)±i^p_,nA^)dP-, (21) 



Mipn,p-) = exp - 2pd)' - 2(pb)^)J 

/ + 00 
M*{pD,p'AM{pD,P~)dpD = 2nqE6{p- - p'„ 
"OO 
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The same relation is valid for functions with (±) indices above. 

The saddle points 7r„ = —2{qExQ — Pd) give the main contribution to the integrals ( pT]) 
with +V'p-,n,T-(a;) functions at xq —>■ ±cx3 (this was first found in Ref. [jT8[ for the d = A case 



without magnetic field). Since the plane wave vanishes in ^ipp_,n,r{x) at 7r_ —>■ ±00, relations 
(pip reduce at Xq ±00 to the space-time uniform field case and can be presented by use 
of formulas from ||l9l as follows 



^IppcnA^) = {'Pf^l^ + ^) +0pi3,n,r(a;) , Xq ±00, (22) 

;0„,,(x°) = CD,_i[±(i-Oe], 

i= (qEx^ - Pd) /\lqE, v = iX/2, C = {A-nqEY^''^ e^^{{-K /2 + i\n2)\/ 



where the function (j)n,r{,x±) was defined in ([TT]), and, Dy{z) is the Weber parabolic cylinder 
function Such solutions were discussed in One can now verify the orthonormality 



relations of the sets{+ippj^^n,r{.x)} and {~''Ppr,,n,r{.x)} . Using transformation (|20| ) one gets the 
orthonormality relations of sets |+V'p_,n,r(a^)} and |~'?/'p_,n,r(a;)| as well. Using the explicit 
forms of solutions (|I6D , (piSD , and relations (0) one can derive the following relations, 

+^(x)G (+1+)^ = +^(x) - Xx)G (+r)^ , 

_7A(x)G (_|-) = -^(x) -+ i;{x)G (+|-) . (23) 
By means of the latter one can get the orthonormality relations for both sets of the solutions, 

[^'4^p-,n,r, ^Vp'_,n',r-') = S(;(;'6rr'Snn'S{p- - p_) , (24) 

where (, (' = ±, 5„„/ is the Kronecker symbol for the discrete spectrum and the 5-function 
for the continuous one. Here r = r' = +1 if = 2, 3. 

To solve a problem of the (±) classification of the solutions one needs to study their 
asymptotic behavior at xq ±oo. From the asymptotic forms (p2D it follows that the 
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asymptotic of the quasienergies of these solution: +e = —qEx*^, is positive and ~e = —qEx^ 
is negative. Thus, the solutions ^%ljp^^n,r{.x) describe particles at Xq — oo, and the solutions 
~i^PD,ri,r{,x) describe antiparticles at xq +00. Since the solutions -ipp_.nT{,x) are orthogonal 
to j^iljp^^n,r{,x), they describe antiparticles at xq —00, and solutions ~^ipp_,n,r{x) describe 
particles at Xq +00 since they are orthogonal to ~il'p-,n,r{x)- One can verify this by taking 
into account that the main contribution to integrals (^) at Xq ±00 for 'til'p,r{x) is given 
by point 7r_ = 0. In this limit the contribution of the plane wave does not depend on x and 
the results of transformation (^) are proportional to a superposition of the solutions in a 



constant uniform field [20 



-0Pi3,n,r(x) = t0„,r-(a;°) exp 0n,r(a;±)i'+ir, 



(25) 



where anr are some coefficients dependent on the plane-wave form. From the asymptotic 
representations (|^) one can see [|l^] that the asymptotic of the quasienergies of these solu- 
tions: = qEx^ is positive and _£ = qEx^ is negative. 

Both sets {±ipp_,n,r{x)} and {'^ipp_^n,r{x)} are orthonormal and complete at any time 
instant. The form of the commutation function, which will be present in the next section, 
can serve as direct proof of the last statement. In c? = 4 such complete sets of the solutions 
were first found in P] . If the plane wave is absent and E = one can get the solutions with 
defined energies from ^ippjy^n,r{x) (|2^ ) considering the next limit in ^^cpn^ri^^)'- 



I =p^nr I 



[^ErrriT^nr + Pd)] ^ GXp [±i l^e^rl x'^) , 

p'jj, d is even 



[(d-l)/2] 



0, d is odd 



where uj is defined in 



In this case one can choose '^4'pjy,n,r{x) =+ ippjy,n,r{x) and 
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Let us select some important properties of the solution. One can see that the matrix ele- 
ments G (cF) diagonal with respect to continuous quantum numbers and spin quantum 
numbers, 

^ Gl^X-,n,r-y_,n',r' = ^rr'S{p- - p'_)S{pi - p[) . . . 5(p2[(d-l)/2]-l " P2[{d-l)/2]-l) 9 (cl^')nn' ' 

(26) 

The solutions ^'4'p_,n,r{x) satisfy the orthonormality conditions on the null-plane, 

J +i>p_,n,ri^)^- +^p'_,n',r'ix)dx+dx^...dx^~^ = 5rr'Snn'5{p^ -p'-), T^T- > 0. (27) 

Since '^i'p_^n,r{x) = for 7r_ < and relations (|^) take place, one gets the following 
representations for G (+|^) G (+|^)^ and G (+|^)^ G (+|^) matrices: 

(G (+1") G = / +ipp_^riA^)'^- +ipp'^,n',r'{x)dx+dx^ ...dx^'^ , 7r_ < 0, 

V ^ ^ ^ ^ / p.,n,r,p'_,n',r' J 

( G (+1^)^ G (+1^) ) = [ ~iljp_,n,r{x)'E.^ ^tljpi_n' r'{x)dx+dx^...dx^^^, 7r_ > 0. (2^ 

V ^ ^ ^ W p_,n,r,p'^,n',r' J 

Thus one can calculate +1? = (? (+|^) (+|^)^ and = (7 (+|^)^ (7 (+|^) matrices using the 
following integrals with the solutions of the squared Dirac equation: 

-Vnn' = J +<Pl_,n,+lA^) ;0p_,n',+l,r(a;)rfx=^C?X^..rfx2[('^-^)/2], ^TT^ < 0. (29) 

III. GREEN FUNCTIONS 

The perturbation theory with respect to the radiative interaction for the matrix elements 
of the processes also has the usual Feynman structure in an external field creating pairs 
|^J2T|J^ . The Feynman diagrams have to be calculated by means of the causal propagator 

S'^{x,x') = c~^i < 0,out\Ttp{x)%jj{x')\0,in >, =< 0, out\0,in > , (30) 

where iplx) is quantum spinor field in the generalized Furry representation, satisfying the 
Dirac equation (|I|), \0,in > and \0,out > are the initial and the final vacuum in this repre- 
sentation, and c„ is the vacuum to vacuum transition amplitude. The propagator S'^{x,x') 
obeys the equation 
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(P^7^ - m) S%x, x') = -^^'^^ {x - x') , (31) 

and is a Green function of the equation. Another important singular function is the com- 
mutation function 

S{x,x')=i ip{x)^%l){x') ^. (32) 

It obeys the homogeneous Dirac equation (|lD and the initial condition 

^(^,^')Lo=.i = ^7°^(x-x'). (33) 

The commutation function S{x, x') is at the same time the propagation function of the Dirac 
equation, i.e. it connects solutions of the equation in two different time instants. 

QED with unstable vacuum has a number of peculiarities. Thus, for instance, in the 
calculation of the expectation values and the total probabilities Green functions of different 
types from (|30D appear ||5|j2l|j3[] : 



S1^{x,x') = i < 0,in\Ti;{x)ip{x')\0,in >, 
S~^{x,x') = i < 0,in\ilj{x)ilj{x')\0,in >, 
S^{x,x') = i < 0,in\ilj{x')ip{x)\0,in >, 
Sf^{x,x') = i < 0,m|^(x)V^(x')T|0,m >, 

S^^^{x,x') = i < 0,out\T^{x)ilj{x')\0,out >, (34) 

where the symbol of the T-product acts on both sides: it orders the field operators to the 
right of its and antiorders them to the left. The functions S^^{x,x'), S^^^{x,x') obey Eq. 
(0), S^{x,x') satisfies Eq. (0) and Sf^{x,x') obeys the equation 



iV,r - m) SUx, x') = 5^''\x - x') . (35) 

As well, all these different kinds of the Green functions are used to represent various matrix 
elements of operators of the current and energy-momentum tensor, and effective action 
beginning with zeroth order with respect to radiative interaction. 
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Solutions (|T6D and (|T8]) with quantum number p_ are especially adapted to calculate 
all the Green functions. One can express the Green functions via solutions (|T6D and (|I8|) 



(x, x') =e{xo- x'o) S- {x, x') - e (x'o - xo) S+ (x, x') , (36) 
S {x, x') = S- (x, x') + S+ {x, x') , (37) 

-OO r / -1 Tll^ 

nr{n . | 



-OO 



S^{x,x')=i dp_ -i^P-,n,rix) g{-\) , ll'p.,n',rix') , (38) 



nr{n'j } 



Sla (a;, x') =e{xo- x'q) (x, x) -0{x'q- Xo) S^^ {x, x') , (39) 
Sin (x, x') = ^ (a^o - Xo) S- {x, x') -e{xo- x'o) S+n {x, x') , (40) 

/ + 00 
dP-Yl ±i^P-,n,r{x)^tpp_^n,r{x'), (41) 
-OO ^j. 

Sout (x, x') =0{xo- x'q) S~^t {x, x') - e (x'o - Xo) Stut {x, x') , (42) 

/+00 _l_ _ ^ 

dp-Yl i^P-,n,r{x) ijp_,n,r{x') . (43) 
-OO „^ 

where all p'j = pj^ the symbol J2nr means the summation over all discrete quantum numbers 
Tij, Tj and the integration over all continuous pj, and the symbol Yl,{nj} means the summation 
over all discrete quantum numbers rij only. Using the relations between the Green functions 
and between the matrices G (cl^ ) "^^^ ^^'^ present the functions S*^, 5*^^ and as follows 

±S^{x,x') = S\x,x')±9{t{xo-x'q))S{x,x') , (44) 

±SUx,x') = SUx.x')±e{T{xo-x',))S{x,x') , (45) 

±SZ,,ix,x') = SUx,x')±eiTixo-x',))Six,x') , (46) 

S^ix, x') = S%x, x') - 5"(x, x') , (47) 

SUx,x') = S%x,x')-S^ix,x') , (48) 
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+00 



dp- J2 -^P-,n,r(x) U+l )gi-\ ) ^1 



(49) 



nr{n'A 



+00 



SP{X,X')=1 dp_ J2 ^^P-,n,r(a;) ^(+1 + ) ^^(+1 ) , ^p_,n'A^') " i^^) 

nr{nj} 

To calculate all the types of Green functions it is sufficient to take sums in S'^{x,x') and 
S°'''''{x,x') only. That will be done below. 
It follows from (|T8|) 

"+CXD 



5" (x,x')= / Q {^-K_)^Y{x,x',p^)dp_. 

-Y{X,X',P_) =iJ2+i^p_n,r{x) +i^p_^n,r{x') ■ 
nr 

Also, taking into account Eqs. (0) one gets 

/ + CO 
e(-7r_)^r(x,x',p_)dp. 
-00 
r+00 

SP{x,x')= e{n_yY{x,x',p_)dp_, 



(51) 
(52) 



(53) 



Due to the fact that +ipp__ n,r {x) and ipp_ n,r {x) have similar form ([T6|), the sums in (^ 
can be taken in the same way. Since 



E 



if (i > 3, and f = S+ if ci = 2, 3, 



the summation over the spin quantum numbers can be done in (E2p to get 



7° + (7±^± + m) 



7r_ 



exp ( -i-a^'T^^ +a 



7° exp 



ln(^7r_) + (In (^vr^ 



/^lP'>{x,x\p.), (54) 

:/(°) (X,X ,p_) = i^;¥?p_,„,r(x) +f*p^,n,rix') , (55) 
n 

+V5p_,n,r (a;) = exp |-^ln(=F7f-)| +0p_,n,r(a;°, X^)0„,r(x_L), 

(56) 



7r 



a = {2qEy^ [{in [T^i)) - In (t^. 



P±M = Oif /U = 0,D, V±^ = V^ if A; = 1,...,D-1 
. 9 



qA^{x'). 
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It is convenient to make the replacement 

in the functions ^K{x_) and 7 J(x_), and the one 



In 



TT 



In r 



in the functions _^^K*{x'_) and _^J*{x'_) (from +0p_,n,r(a; , x ) and 



x'^) defined 



in We do not discuss the integration paths over the variable u, since it is natural to 

consider the plane- wave potential to be an entire function on x_ , in that case the 

integrals ^K{x-), ^J(x_), '^K*{x'_) and ^J*{x'_) are well defined by the integration limits. 
For different forms of the plane-wave potential the integration paths over u are no longer 
arbitrary, however, they easily are extracted from the integration paths shown in Fig.l and 
Fig. 2. Then, one can present some combinations involved in (|55| ) in the following way: 

+K{x.)- = / (;a) + 2 (e''?^+'^ - l) j'^^e'^'^^^qf {x. (u)) du, 

-K + IK* = / (;a) + 2 (e'^^^" + l) e'^^^'g/ (s_ {u)) du, (57) 

-J (x_) - -J* = $ (;a) + 2 Jj'' qf (x_ (m)) e^'^^^tiw gF e'^^^^'g/ (x_ (m')) rfw', 

where 



qf (x_ (m)) [g/ (x_ (m)) + gF/ (m)] rfw. 



l{u) =2 e^^^^^-^'^g/ (x_ («')) rfw', 
Jo 



(58) 



-cx) -iTr/ (2gF)e (±7r: 



1 -exp{-2gFM} 

X-[u) = x_ + y- 



1 — exp |— 2gi?+a| 

The summation over nj in ( |55D can be performed using the Moller formula . After integra- 
tion over pj in (|55|) , applying the operators exp |— i^.ft' (x_) vr^j and exp |— (x_) vr^j , 
and using relations (|57|) , we obtain 



_,_/*-°^ (x, x',p„) = exp |igA — 
1 

~2 

1 



2 - 



ln(=F7r_) + (in (=F7r 



a — 
4 



■y+ }> /i(+a ), 



(59) 



h{-a ) = exp |z$(;a) - ^-(y + /(;a))gF^ coth(gF^ ;a)(y + /(;a)) + lygF^Z(s)| (60) 



where 



(rf-2)/2 / X 

/=i \sm(lHjs)J 

(d-3)/2 / ff \ 

Zid) = c,s''/' n -Tirr ' ^ is odd' (61) 



sin(g_f/'js) 



Q = (47r)-'^/2 - 4)/4} 



and 



A = -£{A^^+A^)dx^. (62) 

Here + is a potential of the constant uniform field F^j^^y, and the integral is taken 
along the line. 

Let us remark that after convenient gauge transformation of electric constant field po- 
tentials the functions ^f^'^^ {x,x',p-) obeys the Klein-Gordon equation, 

(^-2z^ - ^qE + Pi - m^j exp " } = 0" (^3) 

Taking into account the relations 

exp (^—iga'^^FoD '^a^ = cosh (qE '^a^ — ia^^ sinh {^qE 



exp y—i'^cr^^F^^ _^aj = cos (^Q'-f^j + iRj sin (^gi^j _,_aj , (65) 



where 



and matrix Rj is defined by (|15D, one can get relation 

exp {-i^a^'F^^ ;a) 7< ;/(°) (x,x',p_) = 

(tvtx + 7gi^^^ (;«)) exp (-^I^^^F^i ;a) ;/» (x, x',p_) . (66) 
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By using formulas (^5|) , (|66D and Eq. (^) one can transform expression (|5^) to the following 
form, 



:f{'jV + m) ^f{x,x',p. 



(67) 



exp 1-^2^ 



du 



exp < — 



1 r 



ln(=F7r_) + (in (=F7r 



:/(°) (x,x',p_) 



(68) 



In the external field under consideration the real proper time is a function of 7r_ and tt'_ 
because the classical equation of motion has a form nZ^dx- = m~^dS. Thus if ?/_ 7^ one 
can transform the p- integration in Green functions into integration over the Fock-Schwinger 
proper time by making the change of the variable, 



_a . 



(69) 



Then, one gets the following representations for the Green functions: 



5^'"'P(x,x') = (7P + m)A^'"'P(x,x') , 
=fA^(x, x')= / f{x,x',s)ds — Q{±.y-) I 

JTc JTc- 



Ta-ri 



/(x, x', s)ds 



A"(x, x') = / f{x,x',s)ds + Q{y^) / f{x,x',s)ds. 
/SF{x,x') = I f{x,x',s)ds + Q{—y^) / f{x,x',s)ds. 



(70) 
(71) 
(72) 



All the contours of the integrals are shown on Fig. ^. The contours Fc and Fi are placed 
below the singular points on the real axis everywhere outside of the origin, and 



/(x,x',s) 



exp ( -t^a^'-F.^s 



+ (n7)7 / ,.ns-2u,df{{nx{u)))^^^r.HqEs) 



du 



E{ny) 



f^'\x,x',s), 



(73) 



f^^\x,x' ,s) = expjigA} — 



qE 



sinh(gi?s^ 



■ exp{— zm^s + 2$(s) 



-i^iy + lis))qFcothiqFs)iy + lis)) + '-yqFl{s)}Z^,), (74) 
where is defined in (|6T|). The function f^^\x,x', s) has two singular points on the 
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Im s 



r 



1 



o 



Re s 



c 



2qE 



r, 



a 




FIG. 3. Contours of integration Fi, Ts, Fc, F^. 

complex region between the contours r^ — Fi and Fq — Fa. They are situated at the imaginary 
axis: Sq = 0, and qEsi = —in. One can transform the contours Fc — F2 — Fi into F and 
r2 + Ts — Fa into F" (see Fig. ^ with radii tending to zero. Since 

Im s 



O 



Re s 





FIG. 4. Contours of integration F, F' 
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/ fix, x', s)ds = if Vi^y" < , / /(x, x', s)ds = if > 



one can rewrite (|70|),(|7l|) and (^) as follows 



S^{x,x') = {■yV + m)A^{x,x'), 
=fA=^(x,x')= / f{x,x',s)ds - Q{±yo) / f{x,x',s)ds, 
S'^'Plx, x') = (7P + m)A°'P(x, x'), 



A"(x,x')=/ fix,x',s)ds + e{-y^) I f{x,x',s)ds 



(75) 



(76) 



AP{x,x')= f{x,x',s)ds + e{y^) f{x,x',s)ds. 
Jv^ Jrj 

One can check that these expressions are valid for arbitrary x and x'. To prove this one need 

to verify that the functions S^{x,x') and S"''P{x,x'), which are presented via integrals (|75D 

and ([f6D, and by means of representations ( ^I]) and (|53|), are the same solutions of the Dirac 

equation for any x and x'. Thus it is enough, to check first that expressions ( ffSD and (|76D 

obey the Dirac equation for any x and x'. Then, one has to prove that the Cauchy conditions 

for distributions ([75| ) and (|76|) coincide at xq = Xg with ( |5lD and ( p3D respectively. As to 

S'^(x,x'), one can use the proper time representations for 5''^(x,x') and 5'(x,x') functions, 

which follow from (|75|), (|36|) and 



(77) 



S''{x,x') = (7P + m)A^(x,x'), 
A'^(x,x')= / f{x,x',s)ds, 

JVc 



S{x, x) = {-yV + m)A(x, x'), 

A(x, x') = sign(xo — x'o) J f{x,x',s)ds. (78) 
One can see that /(x, x',s) obeys the following equations, 

/(x, x', s) = (P' -m'- l^^'^^M^) fix, x', s) , (79) 

lim /(x, x', s) = i6^'^\x - x') . (80) 

Thus f{x,x',s) is the Fock-Schwinger function P^J2B|] . So, representation ( fTTD of causal 



Green function has the well-known Schwinger form |2^. The concrete representation ([78|) of 
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the commutation function has the same form as the general representation pi]. Similar to 



2^ one can select all singularities of this function and see that it is continuous at xq — Xq. 



In d = A case one can transform this representation to the Fock form [|23[. The analysis of 
space-time singularities in the integrals ( |76D over the contour Tl can be done in a similar 
manner to In the d = 4 the representations ( [75| ) - (^) coincide with the ones found in 



P|. To find the proper time representation we used the special behavior of solutions (|16|). 
But one has the different form of solutions ([T7|) if d is odd, E = and all the imaginary 
eigenvalues Hj of the field tensor are not equal to zero. In this case since solutions {\[7\] 
can represent a special case of the d + 1 general form ([161) , one does not need to calculate 
Green functions of the case independently and representations (ffSD, (|77D and ([78|) are valid 
if one puts E = 0, f{{nx)) = and replaces qE/ sinh{qEs) qH(^d-i)/2/ sin (^qH(^a-i)/2) in 
formulas (|73D and (0) (the functions S"''^{x,x') are equal to zero). 



IV. SOME PHYSICAL APPLICATIONS 

All the information about the processes of particle creation, annihilation, and scatter- 
ing in an external field (without radiative corrections) can be extracted from the matrices 
^ (cl''') ®- These matrices define a canonical transformation between in and out creation 
and annihilation operators in the generalized Furry representation 

a\out) = a\m)G (+|+) + b{in)G (_|+) , 

b{out) = a\in)G (+r) + b{m)G (_r) . (81) 

Here a|^|(m), 6|„|(m), a{„}(m), 6|„}(in) are creation and annihilation operators of in- 
particles and antiparticles respectively and a|„|(o-ut),6|„|(o-ut), a^n}{out), b^n}iout) are ones 
of out-particles and antiparticles, {n} are possible quantum numbers. For example, let us 
calculate the mean numbers of antiparticles created (which are also equal to the numbers 
of pairs created) by the external field from the in-vacuum \0,in > with a given quantum 
number p£,,n,r. By using relations (^TJ) and (|2T|) one finds representation of this quantity: 
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pr),n,T 



--< 0,in\bl^^^^^{out)bp^^nAout)\0,in > 



G 



G 



P£>,n,r,p'jy,n' ,r 



P£),n,r,p'j^,n' ,r 



(82) 



G 



G 



M*ipn,P-)M{p'j,,p.)^. 

p-,n,r,p—,n',r ZTKJJZ/ 



(83) 



where the standard space coordinate volume regularization was used, so that 6{pj — p'j 
Sp^^p'_. Then, by using formulas (^) and (pOf ) one gets 

dp^ 



N.. 



-Vr,r,M*{pD,P-)M{p'^,p_] 



2-KqE 



(84) 



PD=p'r) 



If L/) is the length of the correspondent edge of the space box then the maximum wave 
length of the plane wave is 2Ld- Thus, using the Fourier-series expansion of one gets 



N = 



1 /"^^ ^ 
— / Vnn.dp-. 
vn Jo 



(85) 



One can calculate the quantities ^I^nn', given by Eq. (|29|) , taking into account that the 
operator ttj_ is Hermitian, and 

e'^*^^e-'^''^ = exp l^-^K*qFK -i{K- K*) Tr^j . 

Then all the integrals over in form (^9|) can be expressed in terms of the Laguerre poly- 



nomials |25|. In the special case when n = n' one has 

-Vnn = exp { -ttA - Im -K{xA (qf{P-/qE) + qFRe ~K{x^ 



[d/2]-l ^ [d/2]-l 

J: hA Yl L„^(2/.,), 
i=i J i=i 



hj = -\qH,\{{lm -K2,-i{x.)f + (Im -K2j{x.)f), 7r_ > 0. 



Remember that we are discussing the case in which the constant electric field acts for an 
infinite time. However, one can analyze the problem in finite times T = x^^^^ — a;°„, acting 



similar to |T3[. Then the mean numbers of pairs created by the external field Np^^n,r are the 
same if time T is large enough: ^/qE T ^ 1, y/qET ^ A and qET ^ \pd\- 

Summing over the quantum numbers in (^5]), one can find the total number of pairs cre- 
ated from the vacuum. Using standard regularization with respect to the (c? — l)-dimensional 
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spatial volume V(^d-i) and special regularization with respect to time T of acting of a constant 
electric field where / dp£)Np^^n,r = (lETNp^ n,r, one gets 



N = V( 



{d~i) 



n 



n 



Tm^(3{l) E 



J{d)7^nr^^zj7^— exp i -TT— , , 



E, 



(87) 



where n'^^ is the number density of the created pairs for time T, and the coefficient is 
defined by the next formula special case on / = 1: 

{d-2)/2 

(3{l)= n {qHjCoi\i{l'KHj/E)] , d is even , 

/?(/) = — — n {^^J- coth(/7ri7j/E)} , d is odd . (88) 

Here, i^c = Tn'^/\q\ is the characteristic value of a constant electric field strength. This 
quantity does not depend on the parameters of plane-wave field and is the same as the 



number of pairs created in a constant and uniform field [|T3|]. The corresponding formulas 



for the d = A case were first written in [|26|, and, in fact, can be derived easily from the 



Schwinger formulas [22 



By using proper-time kernel f{x,x', s) ( |73D one can construct the d dimensional form of 
the Schwinger out-in effective action 



out— in 



1 

-tr 
2 



dx 



-1 



f{x, X, s) ds> . 



(89) 



It is not dependent on the parameters of the plane- wave field. And it is not amazing because 
the effective action is a function of the field invariants only which do not depend on the plane 
wave. Taking into account formulas (|65|) one can find the trace in (|89|) representation 



[{d-2)/2] 



p{s) = tr jexp (^-i^a'^'F^^s^ | = 2^^^ cosh(qEs) H cos{qHjS 



(90) 



Then, one calculates the probability for a vacuum to remain a vacuum by using Schwinger 
method 



= exp {-2lmTout-in) = exp {-fiN} 

^ /3(/ + i) r , 



(91) 
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where (3{l) is defined by (|88D. This result coincides with the result from |T^, and in d = 4 



the results coincide with Refs. 122 



Let the operator of the current of the Dirac field operator iIj{x) have the form 



(92) 



and operator of metric energy-momentum tensor (EMT) of the Dirac field operator has the 
form 



T, 



-ican I rpcan 



rpcan 



2 V ' 



(93) 



where T™" is the canonical EMT operator. We are going to discuss the following matrix 
elements with these operators: 



< jf, O,out\j^\0,in > , 

< T^j, 0,out\Tf,^\0,in > 

< J/. >™=< 0,in\j^\0,in > , 

< T^, 0,m|T^,|0,m > , 

< j/. 0, out\jf,\0, out > , 



< T^^ >°"*=< 0, out\Tf,^\0, out > . 



(94) 
(95) 
(96) 
(97) 
(98) 
(99) 



Using the Green functions which were found before, one can present these matrix elements 
in the following form: 



< >'= iqti {-ff,S%x,x)} = iq ti {-f^-f^V^A%x,x')}\^^^, , 

< T,, >^= V4 tr { (7^ {v. + V'l) + 7. {V, + V'l)) S\x, x')} 
= iti{B,,A\x,x')]l^^, , 

< >'"=< ]^, - < J^. , 

< 3, =< 3, - < 3, >" , 



(100) 

(101) 
(102) 
(103) 
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^ rp ^OUt_^ rp _ ^ rp p 



= 1/4 {7^ (p. + p':) + 7. + ^';) } l^V, , 

where the Green functions are given by Eqs. (|76D, <^ and the relat ion 



A (x, x) — A (x, x) 
is used. It is convenient to represent < >"'^ and < T^i, >'*'^ as follows: 



- < J, >"=< J, >(i) + < >(2) , 

- < J. >^=< J. >(^^ - < J. >(^^ , 



where 



< >^^^=iq tr{7^7'^P,A«(x,x')} 
<T^, >«=ztr{5^,A«(x,x')}[^^,, , 
A«(x,x') = -\ I f{x,x',s)ds 



(104) 
(105) 
(106) 
(107) 



(108) 
(109) 
(110) 
(111) 



(112) 
(113) 

, (114) 

2 Jr3+r2+r„ 

and all contributions with derivatives of Q{±y^) functions, which are formally divergent, 
are included in terms < >*^^^ and < T^^ >^'^\ The nature of such divergences is connected 



with infinite time T of acting of a constant electric field and was discussed in |13[ (see also 
below). 

The components < >(i'2) and < T^^^ >(i'2) ^an not be calculated in the framework 
of the perturbation theory with respect to the external background or in the framework of 
the WKB method. Among them only the term < was calculated before in c? = 4 



. Only expression (|101|) for < T^j, >^ has to be regularized and renormalized because 
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of the ultraviolet divergences. Expression ( |10(]| ) for the term < >^ is finite after the 
regularization lifting and equal to zero. The terms < >^^^ and < T^^^ >^^^ are also finite, 
and < >'^^^= 0. The terms < >^'^^ and < T^i, >(^^ have to be regularized with respect 



to time T of acting of a constant electric field |T3| and do not have standard ultraviolet 
divergences. That is consistent with the fact that the ultraviolet divergences have a local 
nature and result (as in the theory without external field) from the leading local terms at 
s — s> +0. The nonzero contributions to the expressions < >^'^^ and < T^,^ >(i'2) g^^e related 
to global features of the theory and indicate the vacuum instability. 

At asymptotic region xq = T/2 — > +00 the densities of current and the EMT of particles 
created are 

" 7^^^ ' ^^^^^ 



according to definitions (|96|) - (^). 

Then, using representations (|102|) - (|107|) , and taking into account terms < >^^'' and 



< T^^ >(2) which are uniform, one gets from (|115|) and ([iT 



f; = 2<j,>('\ (117) 

= 2<T,, >(2) . (118) 



Using special regularization with respect to time T of acting of a constant electric field ||T3 
we can interpret these divergent terms and, correct to first order of ^JqET , obtain 



< J, >^^= 2|g|n-5^, (119) 

< Too >"=< Tdd qETrf'' . (120) 

Other components of < T^^y are of the order of \a{^/qET). 

To study the backreaction of particles created on the electromagnetic field and metrics 
one needs the expressions < >*" and < T^i, >*" . We found the form of < >*"=< 

>(2)= 1/2 < and < T^^j >(2)= 1/2 < T^^ . Let us calculate other terms in 
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expressions ( |104|) . By using formulas (|65|) one can find traces in ( p.01| ) and ( p.07| ) represen- 
tations. One has the next nonzero traces: one presented by expression (ROj) and 



tr|7Vexp = tanh{qEs)p{s) , 

tr 172^72^-1 exp (^-t^a^'-F^.s^ | = - tan(gif,s)p(s) . (121) 



Then, nonzero components of < Tn^ and < Tni, >^^^ 



are 



< T^, >'= / T^,{s)ds , (122) 

< T^, >(i)= -i / r^u{s)ds , (123) 

2 7r3+r2+r„ 



T, 



(s) = 6^(s)p(s)/(°)(x, X, s) if /i = z/ 



&d(s) = -60 (s) = ■ /j^ „ . , &2i(s) = ^-1(3) 



sinh(2gEs) ' sin(2gi/js) 
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